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Q\ • Abstract 

(N 

This is the second paper in a series aiming at searching for infrared pumping 

lines for galactic 1612 MHz OH masers. Our paper I is devoted to the 34.6 jjlvh 

—i. ' absorption lines in ISO SWS spectra towards a large sample of galactic OH/IR 

\& • sources. This paper analyzes the 53.3 //m line in the ISO LWS spectra towards a 

similar sample of OH/IR sources. A search with position radius of 1 arcmin in ISO 

Data Archive (IDA) results in 137 LWS spectra covering 53.3 //m associated with 

47 galactic OH/IR sources and 4 ones associated with megamasers Arp 220 and 

NGC 253. (These two magamasers are included for comparison purpose only.) Ten 

Q-f of these galactic OH/IR sources are found to show and another 5 ones tentatively 

show the 53.3 /im absorption while another 7 sources (our group Ul and U2 sources) 

highly probably do not show this line. The source class is found to be correlated with 



H 



* 



the type of spectral profile: red supergiants (RSGs) and AGB stars tend to show 
strong blue-shifted filling emission in their 53.3 //m absorption line profiles while H II 
regions tend to show a weak red-shifted filling emission in the line profile. GC sources 
and megamasers mainly show symmetrical profile in the line core while megamasers 
tend to show an additional absorption tail on the blue side of the line profile. It is 
argued that the filling emission might be the manifestation of an unresolved half 
emission half absorption profile of the 53.3 /mi doublet which might be produced by 
the transitions among the two levels: 2 n 1 / 2 (J = 3/2) and 2 n]/ 2 (J = 5/2) and their 
closely related levels. The 53.3 to 34.6 /im equivalent width (EW) ratio is close to 
unity for RSGs but much larger than unity for GC sources and megamasers while 
HII regions only show the 53.3 //m line. The pump rate defined as maser to IR 
photon flux ratio is approximately 5% for RSGs. The pump rates of GC sources 
are three order of magnitude smaller. Both the large 53.3 to 34.6 /im EW ratio and 
the small pump rate of the GC OH masers reflect that the two detected 'pumping 
lines' in these sources are actually of interstellar origin. The pump rate of Arp 220 is 
32% — much larger than that of RSGs, which indicates that the contribution of other 
pumping mechanisms to this megamaser is important. A handful of non-detections 
of the 34.6 or 53.3 /jm line or both can be explained partly by the genuinely weakness 
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of the OH masers and partly by some other mechanisms weakening the IR pumping 
lines, such as clumpy OH shell or limb filling emission. 
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1 Introduction 



There are mainly three kinds of pumping mechanisms for stable OH maser: 
by infrared photon, by collision or by chemical reaction(Elitzur, [1992]). Har- 
vey et al. ([1974]) favoured the radiative pumping mechanism for 1612 MHz 
OH maser in the envelope of OH/IR stars because they found that the OH 
maser and IR flux were varying together with almost no phase delay in these 
sources. Elitzur ([1992]) had systematically studied the observational facts and 
the theoretical considerations of the pumping mechanism of the OH masers. In 
circumstances where the radiative pumping is favored, the 34.6 and 53.3 /im 
photons had been thought to be the main pumping photons to excite the 
OH molecules from the ground rotational state 2 n 3 / 2 (J = 3/2) to their up- 
per states 2 n 1 / 2 (J = 5/2) and 2 Y[i/ 2 (J = 3/2) respectively. The exited OH 
molecules then cascade down to invert the F — 1 and F = 2 fine-structure 
levels in the ground state. But astronomers had no chance to probe these ab- 
sorption lines until the Infrared Space Observaotry (ISO - Kessler et al. [1996]) 
opened mid- and far-infrared (IR) observations otherwise impossible behind 
the veiling of the atmosphere. Justtanont et al. ([1996]) and Lutz et al. ([1996]) 
were the first to find in ISO spectra the 34.6 /im absorption line associated with 
OH maser sources. The former found it in NML Cyg, a Red SuperGiant (RSG), 
while the latter found it in Sgr A*, the Galactic Center (GC). The next year, 
Sylvester et al. ([1997]) also profitted from the ISO short wavelength spectrom- 
eter and long wavelength spectrometer (SWS & LWS) observations to find the 
IR pumping absorption line at 34.6 /xm together with the OH cascade emission 
lines at 98.7, 163 and 79iim towards another well known RSG: IRC+10420. 
Thai-Q-Tang et al. ([1998]) modelled the pumping of the OH maser in the 
envelope of this star. Their work confirmed the radiative pumping mechanism 
for the stellar 1612 MHz OH maser. The 34.6 /im absorption in another RSG, 
VY CMa, was also reported by Neufeld et al. ([1999]). These previous works 
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call for systematic check for infrared pumping lines of the stellar OH 1612 MHz 
masers. 

In our previous paper (He et al. [2004], here after Paper I), the 34.6 /im ab- 
sorption is found in only 3 RSGs, 2 GC sources and a megamaser Arp220 
and tentatively in another two objects (one HII region and one Mira) while 
it is expected but not detected in another 16 galactic OH/IR sources (these 
are mainly OH/IR stars, Miras and PNe). Several mechanisms suggested for 
explaining non-detections remain open for further analysis. 

In this paper, the 53.3 /xm spectrum of 102 (out of 137) ISO LWS spectra 
associated with 47 galactic OH/IR sources and 4 ISO LWS spectra associated 
with two megamasers are analyzed. Some new 53.3 /im line detection and non- 
detection sources are found. The 53.3 /im line profiles are classified, discussed 
and correlated to SIMBAD object classes and the detection of 34.6 /im line. In 
Sect. 2, the observational data is described. The detectability of the 53.3 /im 
line is discussed in Sect. 3. The dependence of the OH maser efficiency on the 
ISO IR color is presented in Sect. 4. The profiles are classified in section 5. 
Section 6 compares the results of the 53.3 and 34.6 /im lines while the OH 
pump rate of several sources are derived in Sect. 7. A brief summary is given 
in Sect. 8. 



2 Observation and data processing 



ISO had observed more than 100 OH/IR sources during its 26 months mission, 
but not all of them were observed spectroscopically around 53.3 /im. A search is 
performed in the ISO Data Achive (IDA) for associations with OH/IR sources 
based on a large OH/IR source compilation which combines the catalogue 
provided by B.M. Lewis (including Eder, Lewis and Terzian [1988], Lewis, 
Eder and Terzian [1990], Lewis [1992], Chengalur et al. [1993] and Lewis [1994]) 
with OH/IR catalogues from ATCA telescope (Sevenster et al. [1997a]; [1997b]; 
[2001]) and similar catalogues from many other papers (this unpublished 
large compilation contains 1876 OH/IR sources). Chen et al. ([2001]) cross- 
correlated this large OH/IR compilation with the infrared astronomical satel- 
lite low resolution spectral catalogue (IRAS/LRS) and discussed the LRS 
properties of these sources. An extension of Chen's table of OH/IR sources 
with LRS identification was also given in electronical form in Paper I. In the 
present paper, the same OH/IR source compilation is cross-correlated with the 
ISO/LWS spectra database and 137 ISO/LWS spectra covering the 53.3 /im 
region are found to be associated with 47 galactic OH/IR sources within a po- 
sitional error of 1' . Note that the IRAS position is used to do the identification, 
hence the identification results given here can be slightly different from that 
automatically generated by IDA because the SIMBAD positions of the sources 



used by IDA are not necessarily the same as their IRAS positions. The details 
of the cross-correlation are given in Table 1 in which only the 102 spectra (out 
of 137) for Galactic OH/IR sources and the 4 spectra of megamasers are listed. 
The columns in the table are: (1) the source name; (2) target dedicated time 
(TDT, the unique name of each ISO spectra); (3) ISO LWS spectral resolution 
expressed in minimum resolved wavelength separation in [/im] ; (4) continuum 
flux at 53.3 /im and dispersion (5) as the mean standard deviation of all flux 
data points within the 53 ~ 54 /an range. 

The ISO Spectral Analysis Package (ISAP2.1) x , together with LWS Interac- 
tive Analysis software (LIA10.0), is used to process and analyse all these spec- 
tra. Glitches remained in the Auto Analysis Results (AARs) are removed by 
hand carefully. Small memory effect is neglected by directly averaging together 
the two scans. When large memory effect is present or prominent discrepancy 
exits between the up and down scans, the two scans are averaged separately 
to confirm the detection of the 53.3 /im line feature. 

Here is a brief description of the ISO spectra used in this paper. Some of 
the OH/IR sources (12 sources out of 47) have more than one ISO LWS 
associations. For most of these OH/IR sources, all their associated ISO LWS 
spectra have been processed and considered; the only three exceptions are 
IRAS 17424-2859, IRAS 17441-2822 and IRAS 20255+2712, which have lots of 
repeatedly observed LWS spectra and only part of these ISO LWS spectra 
are processed and presented in Table 1. Therefore, only 102 ISO LWS spectra 
associated with the 47 Galactic OH/IR sources are listed in Table 1. 

For the purpose of comparison, the ISO LWS spectra of two OH megamasers, 
Arp 220 and NGC 253 are also processed. The ISO LWS or SWS spectra have 
been discussed by Skinner et al. ([1997]) and Bradford et al. ([1999]) respec- 
tively. Totally 3 ISO LWS spectra for Arp 220 and 1 spectrum for NGC 253 
are found covering the 53.3 /im region (see Table 1). Hence, there are totally 
49 OH maser sources with 106 associated ISO spectra listed in Table 1. 

ISO LWS observation could be performed in one of the four different possible 
modes: LWS 01, 02, 03 or 04, corresponding to raster scanning of the full 
working wavelength range (43 ~ 196.7 /im) or individual small line ranges, 
Fabry-Perot (FP) high resolution scanning of specified wavelength ranges or 
individual small line ranges respectively. Among the ISO spectra used in this 
paper, there are only eight LWS 02 and three LWS 04 spectra while all others 
are observed in LWS 01 mode. For each observing mode, the actual spectral 
resolution depends on the scan speed. For the 106 LWS spectra used in this 
paper, the spectral resolution around 53.3 /im is found to be about 0.0015 



1 The ISO Spectral Analysis Package (ISAP) is a joint development by the LWS 
and SWS Instrument Teams and Data Centers. Contributing institutes are CESR, 
IAS, IPAC, MPE, RAL and SRON. 



Table 1 

OH/IR sources with their associated ISO LWS spectra covering the 53.3 /iin region. 
Groups and superscripts are explained in the text. Two megamaser galaxies are 
included in group 'A'. 



Name 


TDT 


R[/i] 


*53.3 [Jy: 


v[Jy\ 


(1) 


(2) 


(3) 


(4) 


(5) 


Group 


A 








03507 + 1115 65601707° 


0.07 


427 


2.1 


06053 - 


0622 71002509" 


0.07 


15912 


33.3 


07209 - 


2540 73502338" 


0.07 


1890 


3.2 


17424 - 


2852 32600904" 


0.035 5182 


8.9 


17424 - 


- 2859 49801004" 


0.14 


18284 


145.8 


17430 - 


- 2848 50701028" 


0.071 1847 


19.1 


17441 - 


2822 32601008" 


0.035 6342 


18.9 




32701306* 


0.27 


6517 


129.4 




49400302" 


0.03 


5775 


31.2 


17574 - 


- 2403 09901026*" 


0.07 


15908 


21.7 




12500924* 


0.14 


14975 


54.4 


20255 + 3712 01301603" 


0.07 


9832 


22.6 




12600515* 


0.07 


8425 


38.7 




13400331* 


0.07 


8462 


55.6 




18204615 


0.07 


8837 


147.5 




20306815* 


0.069 8915 


60.1 




21002015* 


0.071 8950 


115.9 




22302315* 


0.07 


8778 


48.7 




22302801" 


0.071 8457 


28.3 




35201327* 


0.14 


8103 


78.2 




49602708* 


0.069 8659 


104.9 




52403008* 


0.07 


8669 


38.1 




53000508* 


0.27 


8589 


195.6 




53102508* 


0.07 


8534 


37.1 



(Table 1 continued) 



(1) (2) (3) (4) (5) 

53803008" 0.07 8557 30.2 

54403605 a 0.07 8603 36.0 

55205608' 0.07 8636 50.5 

55903008 0.07 8600 27.3 

57302606' 0.07 8712 51.6 

57902706 a 0.07 8589 27.1 

58702901" 0.07 8703 28.0 

70601601' 0.07 8897 88.3 

76902601' 0.07 8772 57.2 

77602401' 0.07 8972 49.8 

86000201' 0.071 8935 40.4 

Arp220 27800202 a 0.069 112.8 1.2 

64000801 a 0.069 121.2 2.1 

64000916 a 0.068 121.2 2.3 

NGC 253 24701103" 0.069 964.7 4.4 

NMLCyg 34201304" 0.068 1217 5.7 

55000813™ 0.0015 1919 517.2 

55500942"" 1 0.068 1178 3.3 

Group T 

06319 + 0415 87102811' 0.07 1019 3.5 

16342 - 3814 08402827' 0.071 346.5 4.0 

17431 - 2846 46400917' 0.14 3697 61.0 

18348 - 0526 33000316' 0.071 398 4.6 

34000102 0.14 393.1 4.0 

47201318" 1 0.07 412.8 7.0 

19244 + 1115 31601203' 0.071 795.5 4.1 

36401611™ 0.0015 255.9 613.8 

72400312 0.14 807.5 4.8 

72400414' 0.14 794.2 6.0 



(Table 1 continued) 



(1) (2) (3) (4) (5) 

72400605 0.14 790.7 7.0 

72400702 0.071 794.2 2.5 

Group E 

19039 + 0809 31600514 6 0.035 45.6 1.68 

23412 - 1533 20000913 6 0.071 67 2.74 

Group Ul 

17103 - 3702 28901940 0.07 864.7 10.2 

48202504 0.035 857.2 5.0 

48903704 0.035 846.5 4.4 

49601204 0.035 868.5 5.4 

50301504 0.035 859.7 5.7 

51001804 0.035 865.6 8.6 

67100301 0.07 883.5 7.0 

67801601 0.035 885.7 5.8 

17411 - 3154 46901615 0.07 1382 5.0 

18196 - 1331 34000203 0.071 4476 13.1 

19114 + 0002 31900901 0.072 642.3 4.0 

52500806™ 0.035 660.4 2.7 

52500861 0.29 655.1 2.9 

22176 + 6303 09101821* 0.07 350.8 8.8 

60101805 0.067 11285 28.0 

82301120 0.071 11743 85.8 

82301122 0.07 11771 99.7 

82301123 0.071 11792 39.5 

Group U2 

10197-5750 10300135 0.28 634.7 5.1 

17439 - 2845 69601311 0.27 1854 15.6 



(Table 1 continued) 



(1) (2) (3) (4) (5) 

Group U3 

01037 + 1219 57700513 0.07 163.7 3.5 

57701103 0.07 170.4 2.4 

01304 + 621161300914 0.07 141.8 3.4 

05506 + 2414 83901512 0.07 112.3 2.5 

15452 - 5459 48800916 0.07 339.5 9.3 

17150 - 3224 32702239 0.14 294 5.1 

18050 - 2213 33100802 0.071 290.5 4.4 

18198 - 1249 47801312 0.14 63.9 1.4 

18498 - 0017 32300501 0.071 173.7 4.6 

18560 + 0638 70900321 0.071 128.7 3.3 

19343 + 2926 35501620 m 0.14 86.8 2.9 

36701902 0.071 122.2 4.6 

52000720™ 0.14 131.4 3.7 

52000845 0.071 126.6 2.6 

22177 + 5936 28300920 0.071 94.2 3.0 

Group U4 

07027-7934 14101004 0.14 45.8 2.5 

56700908™ 0.035 29.2 2.3 

56700981 0.28 43.3 2.1 

72200934 0.072 24.4 17.2 

10580 - 1803 22800816™ 0.071 78.7 4.4 

16280 - 4008 08402635 0.071 107.2 4.6 

18272 + 0114 14900719™ 0.035 93.7 3.1 

18437-0302 71501714 0.07196.9 4.1 

18596 + 0315 49901207 0.14 32.5 1.4 

19219 + 0947 54700310 0.07 44.7 2.9 

19255 + 2123 17600528 0.071 48.6 2.4 



(Table 1 continued) 



(1) (2) (3) (4) (5) 

20077 - 0625 34400719 0.072 207.8 9.1 

22036 + 5306 54800798 0.07 119 6.2 

Group U5 

17463 - 3700 32400610 0.27 6.4 4.0 

20000 + 4954 26300417™ 0.14 2.2 1.8 

Notes for column 2: 



The down scan of 09901026 is of bad quality, only its up scan is used. There are data from different 
rasters for 09101821, 63300601, 63300602, the data of raster 3, 1, 1 are used respectively. 

'm': Several spectra arc marked by 'm': 52500806, 47201318, 35501620, 52000720, 56700908, 22800816, 
14900719 and 26300417 were observed in LWS 02 mode; 55000813, 55500942, 36401611 were 
observed in LWS 04 mode. The other spectra were observed in LWS 01 mode. 

'a, t, e': Spectra marked by 'a' means the 53.3 /im line is detected, by 't' means only tentatively detected 
and 'e' means dominated by spurious broad emission feature around 53.3 (im. 

(2 spectra) or 0.068 /im (1 spectrum) for LWS 04 mode, 0.035 (3 spectra), 
0.07 (2 spectra) or 0.14/xm (3 spectra) for LWS 02 mode, 0.035 (10 spectra), 
0.07 (66 spectra), 0.14 (12 spectra) or 0.28 /im (7 spectra) for LWS 01 mode. 
Details of the ISO LWS instrument and observations can be found in Clegg 
et al. ([1996]). 



3 Detectability of the 53.3 /im absorption line 

3.1 Statistics upon the spectra 



Those spectra with the 53.3 /im line feature detected and tentatively detected 
can be defined as group 'a' and 't' and are marked in Table 1 by superscripts 
'a' and 't' respectively; the other non-detections can be defined as group 'u'. 
Another two spectra with spurious broad emission feature around 53.3 /zm 
are defined as group 'e' and marked by superscript 'e' in Table 1. The tenta- 
tive detections are defined as those spectra with measured line depth/height 
smaller than 3 times of the sigma noise level of the continuum flux. The spu- 
rious emission-like feature in group-e spectra is, according to the ISO LWS 
handbook, caused by a near-infrared light leakage in the blocking filters lo- 
cated in front of the detectors. Among the 106 ISO LWS spectra, there are 
21 group-a, 24 group-t, 2 group-e and 59 group-u ones. The group-u can be 



further divided according to their 53.3 /im continuum flux 3 sigma noise level 
into 5 sub-groups: ul, u2, u3, u4 and u5 with the 3 sigma noise level < 2%, 
< 4%, < 10%, < 24% and > 24% of the continuum flux level respectively. 
Totally there are 15 group-ul, 9 group-u2, 16 group-u3, 14 group- u4 and 5 
group-u5 spectra. The group-ul and u2 spectra are considered as highly prob- 
able non-detection cases for the 53.3 /im line while the other group-u3, u4 and 
u5 spectra can be considered as noisy spectra. Therefore, the detection and 
non-detection rates based on the counting of spectra can be estimated to be 
43% and 23% respectively, if the group-t spectra are considered as detections 
and the group-e spectra are excluded. 

Here we do not defined non-detections as we did in paper I where we compared 
the 3 sigma noise level with the expected 34.6 /im line depth estimated from 
the blue peak intensity of the OH 1612 MHz maser because: (1) no theoretical 
hint about the quantitative relationship between the 53.3 /im line and the 
OH 1612 MHz maser line strength is available; (2) the filling emission feature 
widely appearing in the 53.3 /im line profile, as we will see in the following 
discussions, makes it meaningless to compare the directly measured 53.3 /im 
line depth/height with the OH maser line strength. 

The above classification is visualized in Fig. 1 in which the horizontal lines 
delineate how the group-u is further divided into subgroups. This figure also 
shows that the detection of the 53.3 /im line feature generally occurs in spec- 
tra with low 3 sigma noise and high continuum flux, which indicates that the 
selection effect produced by the limited sensitivity of ISO LWS spectrometer 
(sensitivity selection effect) is serious in our sample. Therefore, some weak 
non-detections could be explained by the insufficient sensitivity of the instru- 
ment. Note that two non-detection spectra marked by their source names in 
the figure show abnormally high noise; they are both LWS 04 spectra. Maybe 
the ISO LWS 04 mode is sometimes not as reliable as the other modes. 

Our sample of ISO spectra is biased by selection effects because the obser- 
vations were originally designed for quite different purposes by different ob- 
servers. Beside the sensitivity selection effect mentioned above, another promi- 
nent selection effect in our sample spectra is the different resolutions used for 
different observations. As Fig. 2 shows, most of the spectra were observed 
with a resolution of 0.07 /xm while a non-negligible number of observations 
are performed with higher or lower resolutions (resolution selection effect). 
The different resolutions will affect the strength of spectral features and hence 
affect the detectability of the 53.3 /xm line. 

The detectability of the 53.3 /im OH line feature depends on several facts such 
as noise level, resolution and the genuine strength of the line feature itself. As 
described above, our sample of spectra are affected by instrumental sensitivity 
and resolution selection effects. If the spectra either with 3 sigma noise level 
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Fig. 1. The relative 53.3 /im continuum 3 sigma noise level against the continuum 
flux. The horizontal lines delineate how the group-u spectra are further grouped 



into ul, u2, u3, u4 and u5. 
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Fig. 2. Stack bar figure of spectra of different groups with different resolutions. 

higher than 4% or with resolution other than 0.07 /xm are excluded, we can 
define a relatively more homogeneous sample II with 43 spectra in it. Based 
on sample II another detection and non-detection rate can be re-estimated to 
be 74% (32 spectra out of 43) and 26% (11 spectra out of 43) respectively. We 
do not further constrain the detection and non-detection rate upon the effect 
of the genuine 53.3 /im OH line strength because it is difficult to estimate an 
expected line strength for this OH absorption line due to the reasons described 
earlier. 
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Fig. 3. Stack bar figure of OH/IR sources of different groups with different 53.3 /im 
continuum flux. 

3.2 Statistics upon OH/IR sources 



The OH/IR sources associated with these ISO LWS spectra can be grouped 
in a similar way. Those OH/IR sources with at least one group-a spectrum are 
defined as group 'A'; those with no group-a but at least one group-t spectrum 
are defined as group 'T'; those with neither group-a nor group-t spectra are 
defined as group 'LP and they can be further grouped into Ul, U2, U3, U4 and 
U5 if their associated best quality spectrum belongs to group ul, u2, u3, u4 
and u5 respectively. Another two OH/IR sources with their only one spectrum 
showing spurious broad emission around 53.3 /zm are classified as group 'E'. 
As a result, the 47 galactic OH/IR sources are divided into 10 group- A, 5 
group-T, 5 group-Ul, 2 group-U2, 11 group-U3, 10 group-U4, 2 group-U5 and 
2 group-E sources; the two megamasers, Arp220 and NGC253, both belong 
to group 'A'. 

Another detection and non-detection rate of the 53.3 /zm line can be estimated, 
by counting group-A, T, U galactic OH/IR sources, to be 33% (15 sources out 
of 45) and 16%(7 sources out of 45) respectively. Here again, group-T sources 
are considered as detections and group-E sources are excluded. However these 
rates are also affected by selection effects, as decribed below. 

Beside the instrument related selection effects discussed in Sect. 3.1 that have 
significantly reduced the number of useful spectra in our sample, the dis- 
tribution of our sample sources upon their brightness and IR color can be 
another two potential selection effects, just as we discussed in Paper I. The 
statistics of all our sample sources upon 53.3 /xm continuum flux are plotted 
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in Fig. 3. The source counting distribution in this figure seems to decrease 
from low flux to high flux slightly irregularly, which implies our sample is 
incomplete (brightness selection effect). Another indication from the figure is 
that most group-A, T, Ul and U2 sources are located in high flux groups 
only and group-U3, U4 and U5 sources (noisy sources) only concentrate in 
the two lowest flux groups, which implies that the brightness selection effect, 
together with the instrument related sensitivity selection effect, has largely af- 
fected the detection/non-detection rate estimated above. The other potential 
selection effect, color selection effect, can be checked for in the IRAS two color 
diagram (wherein the IRAS color indices are defined as C2i=log(F 2 5/Fi 2 ), 
C*32=log(i ; 6o/-^25)) in Fig. 6 in which our sample sources distribute in almost 
all areas where OH masers may appear. Therefore the color selection effect is 
not a serious problem in our sample. 



4 ISO Infrared color 



An ISO infrared color Ciso can be defined using the 53.3 /xm continuum 
flux -F53 3 and the 34.6 /im continuum flux F^ 6 from Paper I, i.e., Ciso — 
log(F5 33 /F 3 = 46 ). For those OH/IR sources with multiple ISO LWS observa- 
tions, the _F 5 C 3 3 of most of the reliable spectra listed in Table 1 are used to esti- 
mate a mean value. Totally three spectra: TDT 36401611 of IRAS 19244+1115, 
TDT 35501620 of IRAS 19343+2926 and TDT 09101821 of IRAS 22176+6303 
are rejected when calculating the mean continuum flux because their contin- 
uum fluxes are obviously smaller than that of the other spectra associated with 
the same OH/IR sources. The reason for the too low fluxes may be related 
to the ill performance of the ISO instruments. The resulted mean continuum 
fluxes are listed in column (4) of Table 2. Also listed in this table are columns: 
(1) - source name; (2) - SIMBAD object class; (3) IRAS/LRS classification 
defined by Volk and Cohen ([1989]) and given by Kwok et al. ([1997]); (5) - 
blue peak flux of 1612 MHz OH maser {Fq H uc ) and (6) - reference code for 
OH peak flux and some notes. 

The relation between the blue peak flux of 1612 MHz OH maser {Fq H uc ) and 
Ciso is plotted in Fig. 4 in which Fq H uc is divided by a mean ISO IR flux 
Ff R defined as the mathematic mean of F% 4 6 and _F 5 C 3 3 to remove the effect 
of different distance to the sources. (When F| 4 6 or F 5 C 3 3 is unavailable, the 
other quantity is used for Ff R .) This maser to IR flux ratio can be considered 
as some kind of measure of OH maser 'efficiency'. Most of our sample OH/IR 
sources, irrespect to whether the 53.3 /xm line is detected or not, distribute 
along a uniform inverse correlation between maser efficiency and Ciso except 
five outliers labelled by their source names in the figure. A line is fitted to the 
data not labelled in fig. 4. 
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Table 2 

ISO 53.3 /im and OH 1612 MHz maser spectral quantities of OH/IR sources. These 
objects are grouped in the same way as Table 1. E= 3 3 in this table is a mean value 
for each OH/IR source. 



Name Class 


LRS F 5 C 3. 3 [Jy 


1 rpP,blU( 


: [Jy] ref 


(1) (2) 


(3) 


(4) 


(5) 


(6) 


Group A 










03507 + 1115 Mira 


E 


427 


2.2 


tLVHW 


06053 - 0622 Comp 


H 


15912 


2.03 


tL-HPhD 


07209 - 2540 PV* 


E 


1890 


200 


tLVHW 


17424 - 2852 V* 


I 


5182 


0.12 


tLVHW 


17424 - 2859 IR 


H 


18284 


0.36 


tLVHW 


17430 - 2848 CI* 


H 


1847 


0.95 


LWH92 


17441 - 2822 IR 


- 


6211 


1.39 


ATCAb 


17574-2403 HII 


H 


15442 


34.7 


tLVHW 


20255 + 3712 HII 


H 


8716 


0.22 


AOe(fig) 


,4rp220 IntG 


- 


118.4 


31.4 


Baa87(fig)(lp) 


NGC 253 SeyfG 


P 


964.7 


— - 




NML Cyg V* 


- 


1438 


638 


(priv) 


Group T 










06319 + 0415 IR 


H 


1019 


0.26 


tLVHW(lp) 


16342 - 3814 pAGB 


H 


346.5 


7.1 


SAP93 


17431 - 2846 IR 


H 


3697 


0.21 


tLVHW 


18348 - 0526 Mira 


A 


401.3 


195 


tL-H 


19244 + 1115 pAGB 


E 


796.4* 


45.6 


tLVHW 


Group E 










19039 + 0809 Mira 


E 


45.6 


4 


tLVHW 


23412 - 1533 Mira 


E 


67 


0.06 


ISH94(lp) 


Group Ul 










17103 - 3702 PN 


L 


866.4 


0.2 


ZLP89 


17411 - 3154 OH/IR 


, A 


1382 


224 


tL-H 


18196 - 1331 YSO 


A 


4476 


2 


tL-HPhD 
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(Table 2 continued) 



(1) (2) (3) (4) (5) (6) 

19114 + 0002 PAGB H 652.6 3.98 L89 

22176 + 6303 HII H 11648*0.34 BLS90(fig)(lp) 

Group U2 

10197 - 5750 WR* H 634.7 44 tLVHW 

17439 - 2845 OH/IR H 1854 4.1 tL-H 

Group U3 

01037 + 1219 OH/IR E 167.0 48.2 tLVHW 

01304 + 6211 Mira A 141.8 32 tLVHW 

05506 + 2414 HH H 112.3 0.18 AOa(fig) 

15452 - 5459 pAGB U 339.5 12.3 tL-H 

17150 - 3224 pAGB H 294 2.98 tL-H(lp) 

18050 - 2213 SRPV* E 290.5 16 tLVHW 

18198 - 1249 M1C1 A 63.9 7.4 tLVHW 

18498 - 0017 Mira H 173.7 29.5 tLVHW 

18560 + 0638 Mira A 128.7 20 tLVHW 

19343 + 2926 Em* H 126.7* 0.45 tLVHW 

22177 + 5936 V* A 94.2 41 tLVHW 

Group U4 

07027 - 7934 PN H 35.7 8.62 tL-H 

10580 - 1803 SRPV* E 78.7 2 EN79 

16280 - 4008 PN L 107.2 0.13 LC96(fig) 

18272 + 0114 V*0 - 93.7 0.069 AOe(fig) 

18437 - 0302 IR - 96.9 0.80 BLW94(lp) 

18596 + 0315 OH/IR - 32.5 6.64 ATCA3 

19219 + 0947 PN H 44.7 5.8 tLVHW 

19255 + 2123 PN H 48.6 0.4 tLVHW 

20077 - 0625 V* E 207.8 11.0 tL-H 

22036 + 5306 OH/IR U 119 0.7 tLH91 
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(Table 2 continued) 



(1) 



(2) (3) (4) (5) (6) 



Group U5 

17463 - 3700 PN F 6.4 0.24 ISH94(lp) 
20000 + 4954 Mira E 2.2 3.8 tLVHW 

Column 2: CI* = 'Cluster of Stars', Comp = 'Composite object', Em* = 'Emission Line Star', HH = 
'Herbig-Haro object', IntG = 'Interacting Galaxies', IR = 'Infra-Red source', Mira = 'Variable star 
of Mira Cet type', OH/IR = 'Star with envelope of OH/IR type', pAGB = 'Post-AGB Star', PN 
= 'Planetary Nebula', PV* = 'Pulsating variable Star', SeyfG = 'Seyfert Galaxy', SRPV* = 'Semi- 
regular Pulsating Star', V* = 'Variable Star', V*0 = 'Variable Star of Orion type', WR* = 'Wolf- 
Rayet Star', YSO = 'Young Stellar Object', HII = 'HII (ionized) region', M1C1 = 'Molecular Cloud'. 

Column 4: '*': Following spectra have been excluded due to their abnormal flux value when calculat- 
ing mean continuum flux: TDT36401611 (19244+1115), TDT35501620 (19343+2926), TDT09101821 
(22176+6303). 

Column 6: Reference codes are defined as follow while '(fig)' means the peak flux is measured from the 
published mascr profile figure and '(lp)' means it is a single peak mascr. 

AOe: Lewis( [1994]) LC96: tc Lintel Hekkert & Chapman( [1996]) 

ATCA3: Sevcnstcr ct al.( [2001]) LWH92: Lindqvist et al.( [1992]) 

ATCAb: Sevenster ct al.( [1997b]) SAP93: Silva et al.( [1993]) 

BLS90: Braz ct al.( [1990]) ZLP89: Zijlstra et al.( [1989]) 

BLW94: Blommacrt ct al.( [1994]) tL-H: te Lintel Hekkert et al.( [1991]) 

Baa87: Baan & Hachick( [1987]) tL-HPhD: tc Lintel Hckkcrt's PhD thesis ( [1990]) 

EN79: Engcls( [1979]) tLH91: te Lintel Hekkcrt( [1991]) 

ISH94: Ivison et al.( [1994]) tLVHW: te Lintel Hekkert et al.( [1989]) 

L89: Likkcl( [1989]) (priv): private communication. 
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Fig. 4. Blue peak flux of 1612 MHz OH maser against ISO IR color. Outliers are 
labelled by their source names. 
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Fig. 5. ISO IR color against IRAS color. Sources labelled in Fig. 4 are also labelled 
in this figure by their source names. 

Among the five outliers in Fig. 4, Arp 220 is a megamaser and its too red ISO 
IR color may be due to IR emission from large volumns of cold interstellar 
dust. The too blue color of IRAS 20000+4954 and too red color of the other 
three sources may be due to errors in their ISO continuum flux. This point can 
be checked for by comparing Cjso with IRAS color C32 because the wavelength 
of F£ 4 6 is close to that of F25 and the wavelength of F£ 3 3 is close to that of 
Fqq. The correlation of the two colors is plotted in Fig. 5 for sources with good 
IRAS flux data (with quality factor Q>1). The five sources labelled in Fig. 4 
are also labelled in the figure. As seen from the figure, Cjso of most of our 
sample sources obeys a uniform linear correlation with (723 except four of the 
labelled sources (outliers). A solid line is fitted to those good points (unlabelled 
ones and Arp 220) in the figure. A careful check of ISO positions against IRAS 
and SIMBAD positions shows that, for IRAS 18560+0638 and IRAS 20077- 
0625, the too red ISO IR color is due to a positional error of the 34.6 /xm ISO 
observations which results in a too weak 34.6 /im continuum flux. For IRAS 
17441-2822 (=Sgr B2), its positions in Fig. 4 and 5 seem both not very far from 
the correlations, however we note that the ISO color in Fig. 4 is larger than 
the correlation required while the IRAS color is even larger than the ISO color 
in fig. 5, hence this source, if plotted in pump efficiency-IRAS color relation 
(not shown), will drift further away from the expected correlation. Position 
check shows that the 34.6 /im observation is pointed at a position ~ 35" away 
from the IRAS position. On the other hand the ISO beam size of about 80" 
is larger than its IRAS error ellipse of 23" x 6". Therefore the difference in 
observational position and beam size may account for part of the difference 
between the ISO and IRAS colors. However, the reason for the deviation of 
this source in Fig. 4 is still not clear after our simple consideration because 
this source is located near the complicated crowded galactic center. (Interested 
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readers are referred to Goicoechea and Cernicharo ([2002]) for recent works on 
the modelling of plentiful ISO spectra of this source) For the other source 
under the solid line, IRAS 20000+4954, the ISO position agrees quite well 
with its IRAS and SIMBAD positions, but the ISO spectrum around 53.3 /im 
is quite weak and noisy and hence the ISO 53.3 /zm flux is problematic. In 
a summary, the abnormally too red or too blue ISO IR color of three of the 
outliers in Fig. 4 can be explained by bad flux data while for IRAS 17441-2822, 
a GC source, the situation is still unclear. 

The inverse correlation between Fq H uc and Cjso contradicts with that found 
by other authors such as Chen et al. ([2001]) or Sevenster ([2002]). However 
the increase of maser efficiency with IRAS color found by these authors is 
built upon different sample of OH/IR sources. The correlation found by Chen 
et al. ([2001]) is based on OH/IR sources with IRAS/LRS type-A and E only. 
The correlation found by Sevenster ([2002]) is based on a sample of double 
peak OH maser sources and hence OH masers in star formation regions are 
excluded from their sample. As seen from the IRAS two color diagram for our 
sample (see Fig. 6), it includes a large number of very red cool sources which 
may be HII regions. Hence the inverse correlation between Fq H U£ and Ciso 
might be supported mainly by very cool sources such as HII regions. This 
idea is confirmed by Fig. 7 in which the sources are differentiated by their 
IRAS/LRS type, say, the inverse correlation is mainly supported by H type 
sources that are probably H II regions. The decrease of maser efficiency with 
IR color can be naturally explained by the increase of cool interstellar dust 
content around the cooler HII regions but outside of the maser region, for 
they increase the IR flux and IR color but contribute little to the pumping of 
the OH maser. 



5 The 53.3 yum absorption profiles 



The detected 53.3 /im features in group-a spectra are classified into six groups 
by comparing the spectral profiles by eye: (1) red-shifted absorption with blue- 
shifted filling emission (RaBe); (2) blue-shifted absorption with red-shifted 
filling emission (BaRe); (3) asymmetrical absorption profile of which the blue 
side of the line looks steeper than the red side and so the line center is shifted to 
shorter than 53.3 /im (Asymm); (4) broad profile (BroPf); (5) spurious broad 
emission profile (SpurEm); (6) absorption line towards galactic center sources 
or megamaser galaxies, which is centered precisely at 53.3 /zm (after redshift 
correction for megamasers) and show left-right symmetrical broad profiles 
(GCMega). The spectral profile of NGC253 is peculiar, for it shows a very 
broad absorption tail at the blue side of the line, but it is included in the 
' GCmegd group for the convenience of description. 



o 




Fig. 6. IRAS two color diagram for our sample of OH/IR sources. Only good quality 
IRAS fluxes (with Q>1) are used to make color indices. Regions are from van der 
Veen & Habing ([1988]). 
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Fig. 7. The same as Fig. 4 but sources are differentiated by their IRAS/LRS types. 

In order to visualize the characteristics of the different profile groups, some rep- 
resentative spectra from each group except l SpurErrf and l BroPf ', are shown 
in Fig. 8 while all group-a spectra are plotted in Fig. 9 and Fig. 10. The two 
vertical full lines in Fig. 8 mark the range of the 53.3 /im line, the two dashed 
lines inside them mark the center position of the blue filling emission of 'i?a5e' 
type spectra and that of the red filling emission of the l BaR& type spectra 
while another set of dotted lines mark some nearby features shared among 
different groups. The vertical full lines and dashed lines in Fig. 9 and Fig. 10 
are of similar use except that the two dashed lines are used to mark the center 
of the filling emission and the residual absorption. In Fig. 10, the only one 
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Fig. 8. Representative normalized spectral profiles around 53.3 ^m. Some spectra 
are shifted downwards for clarity. The two solid vertical lines mark the wavelength 
range of the 53.3 /im absorption line. The two dash lines mark the peak position of 
the blue filling emission of 'i?aBe' type spectrum and that of the red filling emission 
of the ""BaR^ type spectrum. Another five dotted lines denote some emission- like or 
absorption-like features shared by different types of spectra. 

l BroPf type spectrum is included among ^BaRe 1 type spectra (the light grey 
one (No. 2) in the figure) because its associated OH/IR source 17574-2403 has 
another spectrum (No. 3) belonging to group L BaRe\ 

Fig. 8 demonstrates that, compared to the symmetrical profile of ^GCMegoU 
type spectra, the blue shifted filling emission in group ^RaBe 1 and the red 
shifted filling emission in group '5ai?e' look quite clear. The latter two groups 
also share emission-like features around 54.5 /zm, as shown by several dotted 
lines in the figure, but their profiles are not completely the same. Generally, the 
line features with filling emission component will become much weaker than 
the ones without filling emission component, and hence the detection of such 
weakened line feature may be quite uncertain due to the limited sensitivity 
and resolution of ISO LWS spectrometer. However, the special case of IRAS 
20255+3712 gives us an important chance to confirm the detection of such a 
weak filling emission feature. This source has a lot of spectra (from No. 4 to 26 
in Fig. 10) and the weak red-shifted filling emission feature is seen in almost all 
of these spectra. If the filling emission feature is unreal, it should not appear 
in so many repeated observations towards the same source. The red-shifted 
filling emission feature appearing in another two independent OH/IR sources 
also confirms its reality. Even more strikingly, the spectral profiles of all 'i?a5e' 
type spectra associated with three different OH/IR sources look quite similar 
to each other not only in the blue-shifted filling emission feature but also in 
the whole wavelength range (from 52 to 55 /iin in Fig. 9). 
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Fig. 9. Normalized spectral profiles around 53.3 /im of i RaBe\ ' Asymm 1 and 
' GCMega 1 type spectra. The spectra are numerated at the right side of the fig- 
ure and those spectra belonging to the same OH/IR source are grouped by two 
short lines. Two vertical full lines mark the edges of the 53.3 /im absorption line 
and two vertical dashed lines mark the center of the filling emission and the resid- 
ual absorption. Spectra are shifted upwards or downwards for clarity. From top to 
bottom they are: (1)65601707, (2)73502338, (3)34201304, (4)55500942, (5)32600904, 
(6)50701028, (7)49801004, (8)32601008, (9)32701306, (10)49400302, (11)27800202, 
(12)64000801, (13)64000916, (14)24701103 
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Fig. 10. The same as Fig. 9 but for ""BaRe! type spectra. Most of these spectra 
belong to 20255+3712. Spectra No. (2) plotted in grey dash dot line belong to group 
'BroPf. From top to bottom they are: (1)71002509, (2)09901026, (3)12500924, 
(4)01301603, (5)12600515, (6)13400331, (7)18204615, (8)20306815, (9)21002015, 
(10)22302315, (11)22302801, (12)35201327, (13)49602708, (14)52403008, 
(15)53000508, (16)53102508, (17)53803008, (18)54403605, (19)55205608, 
(20)57302606, (21)57902706, (22)58702901, (23)70601601, (24)76902601, 
(25)77602401, (26)86000201 
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OH/IR sources can also be classified according to the type of 53.3 /im spectra. 
Generally, different spectra associated with the same OH/IR source should 
belong to the same type. This is true for almost all group-A OH/IR sources 
except IRAS 17574-2403 whose only two spectra belong to group 'BaRe 1 and 
l BroPf respectively. But the 'BroPf type spectrum looks strange if compared 
with the other spectrum of this source and its up and down scans show different 
line profiles. Therefore, this spectrum is unreliable. As a result, the 10 group 
A OH/IR sources together with the two megamasers can be classified into 3 
' RaBe 1 type, 3 'BaRe! type, 2 ' Asymrri! type and 4 'GCMegcH type sources 
(Table 3). The profile of the 5 group T sources are also classified in the same 
way, but their classification is unreliable. 

How is the filling emission produced in the l RaBe' and l BaRe' type sources? 
One may argue that it can arise from some kind of emitting OH clumps with 
moving speed different from the absorbing OH clouds. But the amount of 
red-shift or blue-shift of these filling emission features from the line center 
(53.30 /xm) is measured to be about 0.1 /im which is equivalent to a velocity 
shift of about 560km/s. The amount of red-shift and blue-shift are approx- 
imately the same in all spectra. Such a high velocity shift is impossible to 
explain from the expansion of the envelope around AGB stars, RSGs or H II 
regions. A more reasonable explanation is that the 53.3 /im absorption is a dou- 
blet and one component of the doublet is in full absorption while the other 
is in partial absorption or even in emission. The theoretical separation of the 
two components can be found in ISAP line list database to be about 0.09 /im, 
hence the distance from the two components to 53.30 /im is about 0.45 /im. 
Combined with the ISO instrument profile, the limited spectral resolution, the 
noise of the data and the expansion of the circumstellar envelope, this value 
basically agrees with the red-shift and blue-shift of the filling emission. But it 
is not clear how such asymmetric doublet can be produced in the OH maser 
shell. 

The line profile classification is correlated with the IR object class mainly 
given by SIMBAD and included in Table 3. (The table will be further discussed 
in next section.) As it can be seen from the table, among the three 'i?a5e' 
type sources, IRAS 07209-2540 and NML Cyg are well known red supergiants 
(RSGs) while IRAS 03507+1115 is a well studied nearby low mass evolved 
star. Most possibly this is an AGB star because it also shows SiO (Jewell et 
al., [1991]) and H2O (Bains et al., [2003]) masers, and regular near infrared 
variation with a period of 470 days. Therefore this three sources can be con- 
sidered as a group of evolved stellar OH/IR sources. IRAS 03507+1115 might 
be the first AGB star with its 53.3 /im OH absorption line directly detected 
while its 34.6 /im ISO SWS spectra are too noisy to show OH absorption. 
All three 'BaRe' type sources: 06053-0622, 17574-2403 and 20255+3712 are 
associated with known HII regions. Among the two i Asymm y type sources, 
one is associated to a HII region while the other is a stellar cluster which 
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Fig. 11. Equivalent width (EW) of 53.3 /urn absorption line shown against line center. 
Negative EW means absorption. Positive EW of some RSGs is due to strong filling 
emission. The spectra associated with different IR sources are distributed in separate 
areas. The arrow denotes the position of NGC 253. 



might be accompanied by HII regions too. The l Asymrrf type spectra can 
also be considered as l BaR& type spectra with weaker filling emisson at the 
red side of the line profile. Therefore the ^BaRe 1 and i Asymm y type sources can 
be merged into one group: red-shifted filling emission type and hence the five 
OH/IR sources can be considered as a group of H II regions. Almost all spectra 
associated with HII regions are found to show strong ionized [OIII] emission 
at about 51.8/an, this confirms the HII nature of these IR sources. The ISO 
IR Spectral Energy Distribution (SED) peak position given in Table 4 also 
supports the spectral profile type-object class correlation, i.e., the SED of all 
l RaB& type sources peaks at about 20 /im (as evolved stellar OH/IR sources 
do) while that of many l BaR& or l Asymrrf type sources peaks at about 45 /im 
(rather colder than the former). 

In summary, evolved stellar OH/IR sources have 'i?aSe' type spectra showing 
strong blue-shifted filling emission in the 53.3 /im line profile; HII regions 
have l BaR& or ' Asymm 1 type spectra showing red-shifted filling emission in 
the line profile; galaxy scale IR sources have 'GCMeqa! type spectra mainly 
showing a deep absorption line with symmetrical profile precisely centered at 
53.3 /im. If the filling emission can be explained as filling emission in one of 
the 53.3 /xm doublet components, then one can conclude that evolved stellar 
OH/IR sources tend to show a filling emission in the blue component of the 
53.3 /xm doublet, HII regions tend to show that in the red component while 
no prominent filling emission occurs in galaxy related sources. This difference 
may reflect the different physical conditions in the three kinds of different 
astronomical environments. 
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The above correlation between spectral classification and object type can be 
further confirmed by the investigation of the 53.3 /im absorption line center. 
The integrated line flux, center position and equivalent width of the 53.3 /im 
line (-EW53.3) have been estimated for each group-a and t spectrum, but we 
only give a mean value of the spectral quantities for each OH/IR sources 
in Table 4. The relation between EW53.3 and line center is shown in Fig. 11. 
(The other quantities in the table will be discussed in next sections.) Negative 
.EW53.3 means absorption while positive .EW53.3 means emission. In this figure, 
spectra associated with different classes of objects, H II regions, evolved stellar 
OH/IR sources, GC sources and megamasers, gather in different areas. HII 
regions are located just at the blue side of 53.3 /im due to the red-shifted 
filling emission. Evolved stellar OH/IR sources are located at the red side of 
53.3 /im due to the blue-shifted filling emission. GC sources and megamasers 
are located at just around 53.3 /im. The EW^.z of the magamaser Arp220 
is much larger than that of other sources. The .EW53.3 of GC sources and 
NGC253 is also larger than that of HII regions and evolved stellar OH/IR 
sources. The EW53.3 of some evolved stellar OH/IR sources are positive due 
to the strong filling emission. This figure confirms the correlation between the 
line profile classification and the object class given in the upper discussion. 
The very large .EW53.3 of Arp220 indicates that the OH column density in 
this starburst galaxy is much larger than that in the Galaxy. 



6 Correlation between 34.6 and 53.3 /im absorption lines 



The 34.6 and 53.3 /im photons are generally considered as the main infrared 
pumping agency of radiatively pumped OH 1612 MHz maser, and hence the 
absorption lines at this two wavelengths are expected to appear together in 
one and the same OH/IR source. The OH/IR sources with the line features 
detected or tentatively detected or with good quality spectra at both wave- 
lengths are summarized in Table 3. The columns in the table are: (l)-source 
name; (2)-object class; (3)-IRAS/LRS spectral type; (4)-detection or spectral 
profile type of the 53.3 /im line; (5)-detection of the 34.6 /an line and (6)-the 
peak of spectral energy distribution (SED) estimated by observing the stamp 
plots of ISO SWS and/or LWS spectra by eyes. From this table, it can be seen 
that RSGs, GC sources and megamaser galaxies tend to show both 34.6 and 
53.3 /xm absorptions; HII regions tend to show absorption at 53.3 /im while 
their 34.6 /im lines are not seen due to bad quality data; another two stellar 
sources, including one AGB star and one Wolf-Rayet star, show neither lines. 

The two pumping lines can be further quantitatively compared by plotting 
the equivalent width of the 34.6 /im lines (EW34.6 from paper I) against that 
of the 53.3 /xm lines (-EW53.3 from Table 4) for the five Galactic OH/IR sources 
and the two megamasers Arp 220 and NGC 253 of which both pumping lines 
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Table 3 

Summary of the 34.6 and 53.3 /im absorption features of OH/IR sources and mega- 
masers. Short solid lines are used to separate groups of different object types. Only 
those objects with the detection of the 34.6 or 53.3 ^m line or with good quality 
spectral data at both 34.6 and 53.3 ^m are listed in this table. 



name obj. class 


LRS A53.3 


A34.6 Peak \pm] 


(1) (2) 


(3) 


(4) 


(5) 


(6) 


03507 + 1115 AGB 


E 


RaBe 


N 


<30 


07209 - 2540 RSG 


E 


RaBe 


A 


19 


NML Cyg RSG 





RaBe 


A 


19 


06053-0622 HII 


H 


BaRe 


T 


44.5 


17574-2403 HII 


H 


BaRe 


N 


45 


20255 + 3712 HII 


H 


BaRe 


N 


<45 


17424-2852 HII 


I 


Asymm 





45 


17430 - 2848 cluster 


H 


Asymm 


N 


40 


17431 - 2846 IR 


H 


Asymm(t) 


N 


45 


18348 - 0526 Mira 


A 


Asymm(t) 


U 


27 


19244 + 1115 pAGB 


E 


Asymm(t) 


A 


22 


22176 + 6303 HII 


H 


Ul 


N 


45 


06319 + 0415 IR 


H 


GCmega(t) 





<45 


16342 - 3814 pAGB 


H 


GCmega(t) 


N 


45 


17424-2859 SgrA* 


H 


GCmega 


A 


45 


17441-2822 SgrB2 





GCmega 


A 


80 


Arp 220 megamasei 




GCmega 


A 


50 


NGC 253 megamase] 


• P 


blue— tail 


N 


50 


17411 - 3154 OH/IR 


A 


Ul 


U 


35 


10197 - 5750 WR* 


H 


U2 


U 


28 
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'— ': means no available data 

Column 2: Object types from SIMBAD except several coming from literatures: 

03507+1115 (=IK TAU): We identify the star to be an AGB star from the mascr quantities from 
Jewell et al. ([1991]) and Bains ct al. ([2003]) and some well know facts; 
07209-2540 (= VY CMa): Neufeld ct al. [1999] took it as a RSG; 
17424-2852: HII region cross identified by Chen ct al. ([1995]); 
NML Cyg: Justtanont ct al. [1996] took it as a RSG. 

Column 4: 'RaBe', 'BaRc', 'Asymm', 'GCMega' are the 53.3 /im spectral profile types while 'Ul' and 'U2' 
are the 53.3 (im data quality groups defined in Table 1. 't' in a parenthesis means tentative detection 
and hence the profile classifiction is not sure in this case. Spectra of NGC 253 belong to 'GCMega' 
but show strong and very broad blue absorption tail. 

Column 5: Symbols: 'A' and 'T' mean the detection of the 34.6 /^m absorption is sure or tentative re- 
spectively, 'U' means the expected 34.6 /^m line is not detected while 'N' means the spectra are too 
noisy. 

Column 6: The SED peak position is approximately estimated by observing the stamp plots of ISO SWS 
and/or LWS spectra by eyes. 

are at least tentatively detected (Fig. 12) . As seen from this figure, the 53.3 to 
34.6 /im EW ratio is roughly equal to unity (near to the solid line which shows 
where the two EWs equal to each other) for RSGs but quite larger than unity 
for the four galaxy scale sources: Sgr A* (« 18), SgrB2 (« 8), Arp220 (« 6) 
and NGC 253 (« 25). Both 34.6 and 53.3 /im EWs are much larger for Arp 220 
than for galactic sources and NGC 253. Although H II regions are not plotted 
in Fig. 12 due to the lack of 34.6 /im line data, one can take non-detections 
as detections with very weak line strength (below the 3 sigma noise level of 
the 34.6 /im continuum flux). Therefore their 53.3 to 34.6 /im EW ratio can be 
anticipated to be larger than that of RSGs, although their £W53.3 is similar as 
or only slightly larger than that of the latter. The large 53.3 to 34.6 /im EW 
ratio of the galaxy scale sources may be caused by interstellar OH absorption 
in front of a colder infrared emitting background. The large EW ratio of H II 
regions might be explained in a similar way. The fact that all three RSGs have 
their EW ratio close to unity indicates that equal EW of the two lines might 
be a general characteristic of isolated stellar OH 1612 MHz masers while the 
much larger ratio might be a general characteristic of interstellar OH masers. 
The much larger EW of both lines in Arp 220 may be the result of heavy 
interstellar OH absorption, because Arp 220 is a starburst galaxy and huge 
volumns of the OH molecular clouds can exist in star forming regions. 



7 OH maser pumping 



It has been shown in Sect. 3.1 that the 53.3 /im line is detected or tentatively 
detected in 74% of our good quality LWS spectra and in 33% of the OH/IR 
sources with good quality spectra while it is not detected in 26% of good 
quality spectra and in 16% of OH/IR sources with good quality spectra. It is 
important to discuss whether the detections of the line support the radiative 
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Fig. 12. Mean equivalent width of 53.3 /urn absorption line EW53.3 shown against 
that of 34.6 //m absorption line (.EW34.6). The solid line represents equal EW of the 
two lines. 

pumping mechanism and whether the non-detections undermine it. For the 
non-detections, the most natural explanation is that the OH maser may be 
genuinely very weak and hence the pumping line is too weak to be detected by 
ISO LWS spectrometer. But one can see from Fig. 4 that some non-detections 
(several group Ul, U2 and lots of U3, U4 sources) also have strong maser 
emission, therefore the two OH pumping absorption lines should be strong 
too. Maybe the mechanisms discussed in Paper I, such as clumpy OH shell 
and limb filling emission, are still valid to explain the 53.3 /xm non-detections, 
otherwise these OH 1612 MHz masers should not be radiatively pumped. 

For the detection sources (of group A or T), the spectral profile classification 
results may give some hints on the pumping mechanism. All galactic stellar 
53.3 /im detections turn out to show filling emission in one of their hyper- 
fine component, as discussed in Sect. 5. This may indicate that the asymmet- 
ric doublet profile or filling emission is a common phenomenon for galactic 
OH 1612 MHz maser. Such an asymmetric filling emission might be best ex- 
plained by the selective cascade of photons from energy levels higher than 
2 Ui/2{J = 3/2). We also checked the 34.6 /im absorption line profile of the 
three RSGs and found that the two components of the doublet are basically of 
equal strength in their SWS 07 spectra, therefore such an asymmetric 53.3 /im 
doublet is highly probably produced by transitions among energy levels not 
higher than 2 n!/ 2 (J = 5/2). That is to say, transition from 2 Hi/ 2 (J = 5/2) to 
2 U 1 /2{J = 3/2) in some kind of conditions can produce asymmetry in compo- 
nent strengths in the 53.3 /im doublet but not in the 34.6 /im doublet. The fact 
that the filling emission occurs in different components for RSGs and HII re- 
gions reflects the different physical conditions. If this is true, such asymmetric 
doublet profile may provide a strong proof to identify the radiative pumping 
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of the OH 1612 MHz masers from other pumping mechanisms. Further more, 
this filling emission can be considered also as a potential explanation for those 
non-detections of the line at 34.6 or 53.3 /im or both. 



Combining the 34.6 and 53.3 /im data, one can go deeper into the maser pump- 
ing mechanism. The integrated 53.3 /im line flux (-F533) and equivalent width 
(EW53.3) are estimated by using their best quality spectra and the results 
are given in Table 4 (columns (2) and (3) respectively). Also given in the ta- 
ble are: column (4)-the integrated 34.6 /im line flux (F^g) from Paper I (the 
34.6 /xm data of NGC253 is estimated in this paper from its SWS 02 spec- 
trum TDT 37902123); column (5)-integrated flux (blue peak + red peak) of 
the OH 1612 MHz maser emission {Fq 1 ^) obtained from literature (column (7)) 
and derived OH maser pump rates. The pump rate of radiatively pumped OH 
maser has been defined in Paper I as the integrated OH maser photon flux 
devided by the integrated 34.6 /im IR absorption photon flux. Here we can 
define a better pump rate by deviding the integrated OH maser photon flux 
by the sum of the integrated absorption photon flux of both pumping lines at 
34.6 and 53.3 /im (in column (6)). 



It is easy to see from Table 4 that the pump rate of the three well known 
RSGs: 07209-2540, 19244+1115 and NML Cyg, are approximately the same. 
The pump rate of 19244+1115 given here is a little higher than estimated 
using the photon rate data of Sylvester et al. ([1997]), because our 34.6 /im 
integrated flux is 69% larger than theirs. This difference comes from different 
spectra used. (Readers are referred to paper I for detailed comparison.) The 
nearly similar pump rate of the three RSG OH masers suggests that their 
maser pumping mechanisms may be the same - very possibly pure radiative 
pumping. 



The pump rate calculated for the three galaxy scale sources show extreme 
values. The most prominent is the very small 'pump rate' (~ 10 -5 ) of the two 
GC sources: 17424-2859 (=Sgr A*) and 17441-2822 (=Sgr B2). Such a small 
pump rate supports the previously discussed large interstellar contribution to 
the strong 34.6 and 53.3 /xm absorption lines. The pump rate of the star burst 
galaxy, Arp 220, is 0.32 which is much larger than that of RSGs and even larger 
than the theoretically predicted pump rate of 0.25. Such a large pump rate 
indicates that other pumping processes besides the radiative pumping, such as 
collisional pumping by shock waves, must largely contribute to the observed 
megamaser. Lets take the RSGs as a standard template for radiatively pumped 
maser and designate it a mean pump rate of 0.06, then the contribution of 
radiative pumping to the megamaser of Arp 220 is not higher than 0.06/0.32 ps 
19%. This is only an upper limit because, as discussed above, the two IR 
absorption lines are also due to interstellar OH clouds. 
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Table 4 

Mean line parameters of each group-A and T sources and other quantities for 
OH maser radiative pump rate estimation. Fluxes -F533, ^34*6 are i n un it OI " 
[10~ 19 W/cm 2 ], Fqh is in unit of [10~ 26 W/cm 2 ], equivalent width EW is in unit 
of [10 _4 )i/m]. Fq 1 ^ is expressed by (blue component + red component). Negative 
flux means absorption line. 



name 



(1) 



Fits EW 53 .3 



F™ l R Fqh pumprate notes* 



34.6 



(2) (3) (4) (5) (6) 



(7) 



group A 






03507 + 1115 


1.37 


30.8 


06053 - 0622 


-25.33 


-14.7 


07209 - 2540 


3.96 


19.4 


17424 - 2852 


-50.34 


-92.5 


17424 - 2859 - 


-314.79 


-162.6 


17430 - 2848 


-21.72 


-111.1 


17441 - 2822 


-129.3 


-190.5 


17574 - 2403 


-93.2 


-56.2 


20255 + 3712 


-26.52 


-28.8 


Arp 220 


-5.36 


-446.3 


NGC253 


-23.23 


-227.4 


NML Cyg 


-2.08 


-16.5 


group T 






06319 + 0415 


-0.83 


-7.68 


16342 - 3814 


-5.39 


-143 


17431 - 2846 


-31.6 


-80.6 


18348 - 0526 


-2.82 


-66.2 


19244 + 1115 


-1.12 


-10.5 



-17.6 (970 + 781) 0.079 tL-H(fig) 



-27.7 (1.98 + 1.53) 3.56e - 6 LWH92 



-1.10 (3.15 + 2.85) 1.56e-5 ATCAb(fig) 



-1.00 5667** 0.32 baa87(fig) 

-0.72# 
-11.5 (1130 + 24) 0.041 priv(fig) 



.10 662 0.054 SYL97 



*: The reference codes for OH 1612 MHz maser data are the same as in Table 2 except: SYL97: 
Sylvester ct al. ([1997]) 

(fig) means the integrated 1612 MHz OH maser fluxes are derived by measuring from the figure 
published in the paper. 

**: There is only a single broad maser emission feature for Arp 220. 

#: The 34.6 fan line flux of NGC 253 is estimated by this paper from its SWS 02 spectrum TDT 37902123. 
The 34.6 fim line is only tentatively detected in this spectrum. 
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8 Summary 



Based on the statistical analysis of the 102 ISO LWS spectra (out of 137) 
associated with 47 galacitic OH/IR sources and another 4 spectra of two 
megamasers, we discussed the detectability of the 53.3 /zm OH absorption, 
line profile classification and OH maser pumping. The main conclusions can 
be summarized as follow: 

(1) The 106 ISO LWS spectra associated with 47 galactic OH/IR sources and 
two megamasers are grouped according to the detection of the 53.3 /im 
line and 3 sigma noise level into 7 groups: a, t, ul, u2, u3, u4 and u5. 
The detection and non-detection rate of the line based on the spectra 
counting and free of the instrumental sensitivity and resolution selection 
effect is about 74% (32 spectra out of 43) and 26% (11 spectra out of 43) 
respectively. 

(2) The 49 OH/IR sources are also grouped according to the detection of the 
53.3 /im line in their spectra into 7 groups: A, T, Ul, U2, U3, U4 and U5. 
The detection and non-detection rate of the line based on the galactic 
sources counting is about 33% (15 sources out of 45) and 16% (7 sources 
out of 45) respectively. 

(3) The blue peak flux of the 1612 MHz OH maser is found to decrease with 
increasing ISO infrared color Ciso and this trend is mainly supported by 
HII regions (with IRAS/LRS type 'H'). 

(4) The 53.3 /im line profiles can be classified into 6 types: RaBe, BaRe, 
Asymm, BroPf, SpurEm, GCMega. The filling emission appearing in the 
first three types of spectra is argued to be the manifestation of the asym- 
metric 53.3 /im doublet with one of its find structure components in fully 
absorption while the other partially in absorption or even in emission. 

(5) The OH/IR sources with 53.3 /im detection are also classified according 
to the line profile into 4 classes: RaBe, BaRe, Asymm, GCMega. The 
RaBe type sources are all evolved stellar OH/IR sources (2 RSGs and 

1 AGB star) while the BaRe and Asymm type sources (3 sources and 

2 sources respectively) are usually H II regions. That is to say, evolved 
stellar OH/IR sources tend to show filling emission in the blue doublet 
component of their 53.3 /im line, HII regions tend to show that in the 
red doublet component of the line while galaxy scale sources do not show 
such asymmetry in the line profile. IRAS 03507+1115 is the first AGB 
star directly found to show the 53.3 /xm pumping line. 

(6) GC sources show large 53.3 /zm line equivalent width (.EW53.3) but small 
34.6 //m line equivalent width (-EW 3 4. 6 ); RSGs have small EWu.6 and 
-EW53.3; HII regions also have small EW^.z but do not show 34.6 /xm 
absorption; the megamaser Arp 220 has its EW34.6 and -EW53.3 both very 
large. The large EW53.3 of all four galaxy scale sources and large EWus 
of Arp 220 may be all due to interstellar OH absorption. The -EIU53.3 to 
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EWsas ratio is close to unity for RSGs but quite larger than unity for 
the other three kinds of objects: HII regions, GC sources and megamaser, 
which may also reflect a larger interstellar contribution. 

(7) If the OH 1612 MHz masers in some RSGs and H II regions of our sample 
are really radiatively pumped, the filling emission in one of the 53.3 /zm 
doublet components is probably produced by the transitions between 
the two OH energy levels: 2 U 1 / 2 {J = 3/2) and 2 Yl 1 / 2 (J = 5/2) or their 
closely related levels. The fact that the filling emission occurs in different 
components for the two kinds of objects may reflect different physical 
conditions. 

(8) The pump rate of RSGs is about 0.06 while that of GC sources is very 
small and that of Arp 220 is very large. Perhaps the pump rate of RSGs 
represents a typical value for purely radiatively pumped OH 1612 MHz 
masers. Low pump rate of GC sources is due to interstellar OH aborption 
while very large pump rate of Arp 220 indicates a large contribution of 
other pumping mechanisms (e.g., collisional pumping). 

(9) A handful of OH/IR sources with both pumping lines not detected in their 
better quality spectra still need to be explained. The genuinely weakness 
of the OH maser and hence of the two IR pumping lines can explain part 
of these non-detections. Maybe the alternative explanations presented in 
Paper I, such as clumpy OH shell and limb filling emission, are also valid. 



Acknowledgements 



Guest User, Canadian Astronomy Data Center, which is operated by the 
Dominion Astrophysical Observatory for the National Resarch Council of 
Canada's Herzberg Institute of Astrophysics. We thank Dr. R. Szczerba for 
the help in processing ISO spectra and those astronomers who have prepared 
the ISO projects to observe the OH/IR sources used in this paper. Our work 
is supported by the Chinese National Science Foundation under Grant No. 
10073018, the Chinese Academy of Sciences Foundation under Grant KJCX2- 
SW-T06 and the Yunnan Natural Science Fund (2002A0021Q). 



References 

[1987] Baan W.A., Hachick, A.D., 1987, ApJ 318, 139 (baa87) 

[2003] Bains I., Cohen, R. J., Louridas, A., et al., 2003, MNRAS 342, 8 

[1994] Blommaert J.A.D.L., Van Langevelde H.J., Michiels W.F.R, 1994, A&A 287, 
479 (BLW94) 



32 



1999 
1990 
1995 
2001 
1993 
1996 
1988 
1992 
1979 
2002 
1974 
2004 
1994 
1991 
1996 
1996 
1997 
1992 
1994 
1990 
1989 
1992 
1996 
1999 



Bradford CM., Stacey G.J., Fischer J.A., et al., 1999, ESASP 427, 861 

Braz M.A., Lepine J.R.D., Sivagnanam P., et al., 1990, A&A 236, 479 (BLS90) 

Chen P.S., Gao H., Xiong G.Z., 1995, ApJS 100, 389 

Chen P.S., Szczerba R., Kwok S., et al., 2001, A&A 368, 1006 

Chengalur J.N., Lewis B.M., Eder J., et al., 1993, ApJS 89, 189 

Clegg P.E., Ade P.A.R., Armand C, et al., 1996, A&A 315, L38 

Eder J., Lewis B.M., Terzian Y., 1988, ApJS 66, 183 

Elitzur M., 1992, Astronomical Masers (book), Kluwer 

Engels D., 1979, A&AS 36, 337 (EN79) 

Goicoechea J.R., Cernicharo J., 2002, ApJ 576, L77 

Harvey P.M., Bechis K.P., Wilson W.J., et al., 1974, ApJS 27, 331 

He J.H., Szczerba R., Chen P.S., et al., 2004 (submitted to A&A, paper I) 

Ivison R.J., Seaquist E.R., Hall P.J., 1994, MNRAS 269, 218 (ISH94) 

Jewell, P. R., Snyder, L. E., Walmsley, C. M., et al., 1991, A&A 242, 211 

Justtanont K., de Jong Th., Helmich F.P., et al., 1996, A&A 315, L217 

Kessler M.F., Steinz J.A., Anderegg M.E., et al., 1996, A&A 315, L27 

Kwok S., Volk K., Bidelman W.P., 1997, ApJS 112, 557 

Lewis B.M., 1992, ApJ 396, 251 

Lewis B.M., 1994, ApJS 93, 549 (AOe) 

Lewis B.M., Eder J., Terzian Y., 1990, ApJ 362, 634 

Likkel L., 1989, ApJ 344, 350 (L89) 

Lindqvist M., Winnberg A., Habing H.J., et al., 1992, A&AS 92, 43 (LWH92) 

Lutz D., Feuchtgruber H., Genzel R., et al., 1996, A&A 315, L269 

Neufeld D.A, Feuchtgruber H., Harwit M., et al., 1999, ApJ 517, L47 



1997a] Sevenster M.N., Chapman J.M., Habing H.J., et al., 1997a, A&AS 122,79 

(ATCAa) 

1997b] Sevenster M.N., Chapman J.M., Habing H.J., et al., 1997b, A&AS 124,509 

2001] Sevenster M.N., van Langevelde H.J., Moody R.A., et al., 2001, A&A 366, 
481 (ATCA3) 

[2002] Sevenster M.N., 2002, AJ 123, 2772 



33 



[1993] Silva A.M., Azcarate I.N., Poppel W.G.L., et al., 1993, A&A 275, 510 
(SAP93) 

[1997] Skinner C.J., Smith H.A., Sturm E., et al., 1997, Nature 386, 472 

[1997] Sylvester R.J., Barlow M.J., Nguyen-Q-Rieu, et al., 1997, MNRAS 291, L42 

[1989] te Lintel Hekkert P., Versteege-Hensel, H. A., Habing, H. J., et al.,1989, 
A&AS 78, 399 (tLVHW) 

[1990] te Lintel Hekkert P., 1990, PhD. thesis, University of Leiden (tL-HPhD) 

[1991] te Lintel Hekkert P., 1991a, A&A 248, 209 (tLH91) 

[1991] te Lintel Hekkert P., Caswell J.L., Habing H.J., et al., 1991b, A&AS 90, 327 
(tL-H) 

[1996] te Lintel Hekkert P. and Chapman J.M., 1996, A&AS 119, 459 (LC96) 

[1998] Thai-Q-Tang, Dinh-V-Trung, Nguyen-Q-Rieu, et al., 1998, A&A 331, 317 

[1988] van der Veen W.E.C.J., & Habing H.J., 1988, A&A 194, 125 

[1989] Volk K., Cohen M., 1989, AJ 98, 931 

[1989] Zijlstra A.A., te Lintel Hekkert P., Pottasch S.R., et al., 1989, A&A 217, 157 
(ZLP89) 



34 



